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In the framework of quantum computational tensor network, which is a general framework of 
measurement-based quantum computation, the resource many-body state is represented in a 
tensor-network form (or a matrix-product form), and universal quantum computation is performed in a 
virtual linear space, which is called a correlation space, where tensors live. Since any unitary operation, state 
preparation, and the projection measurement in the computational basis can be simulated in a correlation 
space, it is natural to expect that fault-tolerant quantum circuits can also be simulated in a correlation space. 
However, we point out that not all physical errors on physical qudits appear as linear completely-positive 
trace-preserving errors in a correlation space. Since the theories of fault-tolerant quantum circuits known so 
far assume such noises, this means that the simulation of fault- tolerant quantum circuits in a correlation 
space is not so straightforward for general resource states. 

Quantum many-body states, which have long been central research objects in condensed matter physics, 
statistical physics, and quantum chemistry, are now attracting the renewed interest in quantum informa- 
tion science as fundamental resources for quantum information processing. One of the most celebrated 
examples is one-way quantum computation 1-3 . Once the highly- entangled many-body state which is called the 
cluster state is prepared, universal quantum computation is possible with adaptive local measurements on each 
qubit. Recently, the concept of quantum computational tensor network (QCTN) 4 " 6 , which is the general frame- 
work of measurement-based quantum computation on quantum many-body states, was proposed. This novel 
framework has enabled us to understand how general measurement-based quantum computation is performed 
on many other resource states beyond the cluster state. The most innovative feature of QCTN is that the resource 
state is represented in a tensor-network form (or a matrix-product form) 7 " 9 , and universal quantum computation 
is performed in the virtual linear space where tensors live. For example, let us consider the one-dimensional open- 
boundary chain of N qudits in the matrix-product form 



d-i 



(i) 



fci=0 



k N = 0 



where {|0), \d — 1)} is a certain basis in the d- dimensional Hilbert space (2 < d < oo), |L) and \R) are D- 
dimensional complex vectors, and {A[0], A[d — 1]} are D X D complex matrices. Let us also define the 
projection measurement Mq$ on a single physical qudit by 



M B j = {\«ej), |/W),|2>, 1», 



where 



,)scos-|0> + e*sin-|l), 



|A 



04/ 



0. V 

sin-|0>- 



-e t<p cos- II), 



0 < 6 < 7i, and 0 < </> < In. If we do the measurement Mq^ on the first physical qudit of Eq. (1) and if the first 
physical qudit is projected onto, for example, |a^) as a result of this measurement, the state Eq. (1) becomes 
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V (L\A[k N }...A[k 2 ] r£ L \R)\k N , ...,* 2 >®|a^), 
where 

0 0 
A[<x e j\ = cos-A[0] + e~^ sin-A[l]. 

Now we can imagine a virtual linear space where A's, \R), and \L) live. 
This virtual linear space is called the correlation space 4 " 610 " 12 . Then, 
the above equation can be interpreted that the "gate operation" 

\R)-> -n — [ 6, ^\ j j \R) is implemented on the "quantum state" \R) in 

ll A L a ^JII 

the correlation space. In particular, if A[0], A[l], 0, and <j> are appro- 
priately chosen in such a way that A [a^] is proportional to a unitary, 
we can "simulate" the unitary evolution T . — [ 9, ^\ , , \R) of the vector 

\R) in the correlation space. The core of QCTN is this "virtual 
quantum computation" in the correlation space. If the correlation 
space has a sufficient structure and if A's, \L), and \R) are appropri- 
ately chosen, we can "simulate" universal quantum circuit in the 
correlation space 4 " 6 ' 10 " 12 . 

For the realization of a scalable quantum computer, a theory of 
fault-tolerant (FT) quantum computation 13 " 17 is necessary. In fact, 
several researches have been performed on FT quantum computa- 
tion in the one-way model 318 " 22 . However, there has been no result 
about a theory of FT quantum computation on general QCTN. (See 
the recent progress to this direction 23 " 26 .) In particular, there is severe 
lack of knowledge about FT quantum computation on resource states 
with d > 3. It is necessary to consider resource states with d > 3 if we 
want to enjoy the cooling preparation of a resource state and the 
energy-gap protection of measurement-based quantum computa- 
tion with a physically natural Hamiltonian, since no genuinely 
entangled qubit state can be the unique ground state of a two -body 
frustration- free Hamiltonian 27 . 

One straightforward way of implementing FT quantum computa- 
tion on QCTN is to encode physical qudits with a quantum error 
correcting code: 

^■••E (L\A[k N ]...A[k l ]\R)\k N , 

ki=0 k N = 0 

where \kA (i = 1, N) is the encoded version of \k t ) (such as 
1 6) = 1 000; and |l) = |111), etc.) In fact, this strategy was taken in 
Refs. 21 ' 22 for the one-way model (d = 2), and it was shown there that a 
FT construction of the encoded cluster state is possible. For d > 3, 
however, such a strategy is difficult, since theories of quantum error 
correcting codes and FT preparations of the encoded resource state 
are less developed for d > 3. Furthermore, if we encode physical 
qudits with a quantum error correcting code, the parent Hamil- 
tonian should no longer be two -body interacting one. 

The other way of implementing FT quantum computation on 
QCTN is to simulate FT quantum circuits in the correlation space. 
Since any unitary operation, state preparation, and the projective 
measurement in the computational basis can be simulated in a cor- 
relation space (for a precise discussion about the possibility of the 
measurement, see Ref. 12 ), it is natural to expect that FT quantum 
circuits can also be simulated in a correlation space. An advantage 
of this strategy is that theories of FT quantum circuits for qubit 
systems are well developed 13 " 17 . In fact, this strategy was taken in 
Refs. 1819 for the one-way model (d = 2). They introduced a method 
(which we call "the ensemble method" since the ensemble of all 
measurement results are considered) of simulating quantum circuits 
in the correlation space of the cluster state, and showed that all 
physical errors on physical qubits can be linear completely-positive 
trace-preserving (CPTP) maps in the correlation space of the cluster 



state. (For details, see Sec. I of the supplementary material.) This 
means that FT quantum circuits can be simulated in the correlation 
space of the cluster state. 

In this paper, however, we point out that it is not so straightfor- 
ward to simulate FT quantum circuits in the correlation space of a 
general resource state. In order to see it, we consider two standard 
methods of simulating quantum circuits in the correlation space. 
First, we directly apply the ensemble method of Refs. 1819 to the 
one- dimensional Affleck- Kennedy-Lieb-Tasaki (AKLT) state 28 ' 29 , 
as an example, and see that not all physical errors on physical qutrits 
can be linear CPTP maps in the correlation space of the AKLT state. 
Since all theories of FT quantum circuits known so far assume such 
noises 13 " 17 , this means that it is not so straightforward to apply these 
FT theories to quantum circuits simulated in the correlation space of 
general resource states. (Note that quantum error correction is not 
equivalent to fault- tolerant quantum computation: non-CPTP errors 
can be corrected by quantum error-correcting codes, but that does 
not necessarily mean the possibility of fault-tolerant quantum com- 
putation.) We give some intuitive explanations of reasons why the 
cluster state is so special and why not all resource states work as the 
cluster state. Second, we consider another natural method of simu- 
lating quantum circuits in the correlation space, which we call "the 
trajectory method" (since only a specific history of measurement 
results, i.e., a specific trajectory, is considered). This method is 
another standard way of simulating quantum circuits in the correla- 
tion space when the system is assumed to be error-free 4 " 610,28 . We 
show a general theorem which says that if d > 3 not all physical errors 
on physical qudits can be linear CPTP maps in the correlation space 
if we consider the trajectory method. This means that the trajectory 
method does neither work in general. 

In short, we show in this paper that it is not so straightforward to 
simulate FT quantum circuits in the correlation space of a general 
resource state. Since all errors behave nicely in the correlation space 
of the cluster state 1819 , less attention has been paid to the difference 
between a real physical space and a correlation space of a general 
resource state. Our results here suggest that these two spaces can be 
different, and because of the difference, simulations of FT quantum 
circuits can be difficult in a general correlation space. To our know- 
ledge, there is so far no other methods than those two methods, i.e., 
the ensemble method and the trajectory method. However, of course 
there is a possibility that a new method will be developed in a future. 
Therefore, what we show here is not the impossibility of making 
QCTN fault-tolerant. Such a new method might be able to make 
all QCTN fault-tolerant. We hope that our results here will help to 
tackle such a challenging subject of a future study. 

Results 

Assumption. Throughout this paper, we make the following 
assumptions: Since the MPS |¥(L, R))i> Eq. (1), is a resource state 
for measurement-based quantum computation, we can assume 
without loss of generality that A[a^], A[/^], A [2], A [3], ...,A[d 
— 1] are unitary up to constants: 

A[a^] = c a l/ a , 

A[p e j]=cpUp, 

A[2]=c 2 U 2 , 

A[3]=c 3 U 3 , 

A[d-i\=c d -iUd-i, 

where c a , cp, c 2 , ... Q-i are rea l positive numbers, C7 a , Up, U 2 , 
U d - 1 are unitary operators, and 
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A[p u ] = sm-A[0]-e- 



'cos-A[l] 



This means that any operation implemented in the correlation space 
by the measurement Me,<j> on a single physical qudit of |^(L, R)) 1 is 
unitary. Note that this assumption is reasonable, since otherwise 
|*P(L, R))i does not seem to be useful as a resource for 
measurement-based quantum computation. In fact, all known 
resource states so far including the cluster state and the AKLT 
state, satisfy this assumption by appropriately rotating each local 
physical basis. Furthermore, we can take c a , cp, c 2 , . . ., c d - x such that 

d-i 



/f m |DU E (L\A[k N ]...A[k l ]\R)\k N ,...,k 1 
VMI L /> \ R ))k u ...,k N 



h,...,k N 



and we can redefine A[ki\/\/C->A[kj\. 

Simulation on the AKLT state with the ensemble method. Let us 

first consider the direct application of the ensemble method of 
Refs. 18 ' 19 , which was used for the cluster state, to the AKLT state. 
As is shown in Refs. 1819 , all physical errors can be linear CPTP 
maps in the correlation space of the cluster state. (This is also the 
case for the tricluster state 30 . See Sec. Ill and IV of the supplementary 
material.) However, here we show that if we directly apply the 
ensemble method to the AKLT state, not all physical errors can be 
linear CPTP maps in the correlation space of the AKLT state. 

Roughly speaking, in the ensemble method, we consider the mix- 
ture of all measurement results (measurement histories) during the 
feed-forwarding stage. For the cluster state, it works well, since the 
desired gate operation can be implemented deterministically because 
of the feed-forwarding, and if we trace out states which registered 
each measurement results after the feed-forwarding, the entire state 
again becomes pure. (For details, see Ref. 19 and the Supplementary 
Material.) 



The one-dimensional AKLT state is the matrix-product state 
defined by d = 3,A[0] = X,A[1] =XZ,andA[2] = Z. Let us assume 
that a CPTP error p-> ^2J = iFjpFj > where ^2J= 1 FjFj=I, occurs on 
a physical qutrit of the AKLT state. If we simulate quantum circuits 
in the correlation space of the AKLT state after this error with the 
ensemble method, the map 

\R){R\^S AKLT (\R}{R\)^^ S i\ R )( R \ S i 



is implemented in the correlation space, where H/s are certain opera- 
tors which depend on {Fj}. (For details, see Sec. II of the supplement- 
ary material.) For example, let us consider the error with 

w = 1 and F 1 = ^(|«^)+|^))(0|-^(|a^)-|^)) 



(1| + 1 2) (2\. Then, we can show that 



X#a,=«+ai><ii 



where 3 is a certain positive number. This means that the map Saklt 
is not linear CPTP. 

Intuitive explanations. We have seen that the ensemble method of 
Refs. 1819 for the cluster state cannot be directly applied to other 
resource states, such as the one-dimensional AKLT state. Why the 
cluster state is so special? And why it does not work for other resource 
states? Although the complete answer to these questions is beyond 
the scope of the present paper, since the study of QCTN itself has not 
been fully developed (for example, no one knows the necessary and 
sufficient condition for tensor- network states to be universal 
resource states), let us try to give some intuitive explanations here. 

Figure 1 illustrates the reason why all physical errors become 
CPTP maps in the correlation space of the one-dimensional cluster 
state. Let us first consider the ideal case (a) where there is no error. 
The actual protocol (a-1) is mathematically equivalent to the "input- 
output" picture (a-2) where the physical input state | is teleported 
into the left-edge of the short-length cluster chain (indicated in yel- 
low) and finally the physical output state is extracted from the 
right-edge of the short-length cluster chain (indicated in yellow). 
Since both and are physical states, what is going on in the 
correlation space which maps the input state | to the output state 



(a) Ideal case 
(1) actual protocol 

oo ooo- 

input output 

measurement 

(including feedforward) 

oo 



(2) "input-output" picture 

oo 00*0 



\*) ? 



— input output 
measurement 
mathematically ( i nc | udin g feedforward) 
equ-valent ^ ^ 



(b) Imperfect case 
(1) actual protocol 

DO OOO 

input output 
imperfect measurement 
(including feedforward) . 

x bo- 




imperfection 



oo 



mathematically 
... equivalent 



oo 



(2) "input-output" picture 

"do oo ¥ o 

input output 
imperfect measurement 
(including feedforward) 



Figure 1 | The "input- output" picture for the one-dimensional cluster state. 
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can be described by a linear CPTP operation. In other words, if 
we can describe measurement-based quantum computation with this 
"input-output" picture 19 , the map implemented in the correlation 
space is guaranteed to be a linear CPTP operation 19 . 

For the cluster state, this "input- output" picture also holds even if 
there is an error 19 : In the imperfect case, Fig. 1 (b), let us assume that 
the input state is degraded by an error and becomes a mixed state p. 
We also assume that measurements are imperfect. However, we can 
still consider a similar "input- output" picture (b-2) 19 , which corre- 
sponds to the actual protocol (b-1), and again the physical state p is 
mapped into another physical state p\ which means that what is 
going on in the correlation space which maps p to p' can be described 
by a linear CPTP operation. 

Note that two special properties of the cluster state enable such an 
"input-output" picture. First, the one-dimensional cluster state can 
be decomposed into small pieces of one-dimensional cluster states by 
applying nearest-neighbour two-body unitary operations (i.e., CZ 
gates). As we can see from Fig. 1 (a-2) and (b-2), this property is 
necessary for allowing the "input- output" picture. Second, the num- 
ber of qubits that are measured in order to implement a specific gate 
does not depend on the measurement results. In other words, for the 
cluster state, a specific gate can be implemented up to Pauli bypro- 
ducts at a fixed site irrespective of measurement results. Such a deter- 
ministic implementation of a gate at a fixed site is necessary for the 
deterministic (i.e., trace-preserving) "output" in the "input- output" 
picture, since in the ensemble method, the ensemble (mixture) of all 
measurement results are considered: If the site where the desired gate 
operation is completed depends on the measurement results, we 
cannot extract the same output state at a fixed site irrespective of 
measurement results as is shown in Fig. 1 (a-2) and (b-2). 

On the other hand, such an "input-output" picture seems to be 
impossible for the one-dimensional AKLT state, because of the follow- 
ing two reasons: First, as is shown in Fig. 2 left, no nearest-neighbour 
two-body unitary operation can decompose the one-dimensional 
AKLT chain into two chains due to the existence of the non-vanishing 
two-point correlation in the AKLT state. (If the one-dimensional 
AKLT chain can be decomposed into two chains by such a unitary, 
it contradicts to the well-known fact that the two-point correlation is 
non-vanishing in the AKLT state.) Second, we cannot deterministically 
implement a specific gate at a fixed site of the AKLT chain irrespective 
of measurement results 4,5,28 . In short, the "input-output" picture seems 
to be impossible for the AKLT state. If we can no longer use the "input- 
output" picture, it is not unreasonable that we have some anomalous 
maps in the correlation space since the correlation space is not a 
physical space but an abstract mathematical space. 

Simulating quantum circuits with the trajectory method. As we 

have seen, it is not always possible for general resource states to 
implement a specific gate at a fixed site irrespective of measure- 
ment results. This fact prohibits a general resource state from allow- 
ing the "input-output" picture. One might think that if we abandon 
such a deterministic "output" at a fixed site for all measurement 



results, and if we just consider a specific history (trajectory) of 
measurement results, we might be able to avoid the emergence of 
non-CPTP errors, such as Saklt- (Physically, this means that we 
project the system onto a pure state at every measurement step: 

P|<F) 



\\pm 



with the corresponding projector P. Note that considering a specific 
trajectory of certain measurement results does not mean that if we 
obtain other trajectories, such events are discarded as done in, for 
example, linear optical quantum computation. Rather, it means that 
we describe measurement-based quantum computation with respect 
to each trajectory instead of the ensemble of all trajectories.) If the 
system is assumed to be error-free, this "trajectory method" is 
another standard way of simulating quantum circuits in the 
correlation space 4 " 6,10,28 . (Note that in this trajectory method, 
correct unitary operators can be implemented in the correlation 
space if there is no error, although what we are physically doing 
are projections, i.e., non-trace-preserving operations). 

However, we here show that such a natural another way of simu- 
lating quantum circuits does neither work if d > 3. In other words, we 
can show the following theorem. (For a proof, see Methods.) 

Theorem: If d > 3, there exists a single-qudit CPTP error 6 which 
has the following property: assume that 8 is applied on a single 
physical qudit of Eq. (1). If the measurement Mo, 4 is performed 
on that affected qudit, and if the entire system is projected onto a 
pure state as a result, a non-TP operation is implemented in the 
correlation space. ■ 

Discussion 

In order to gain an intuitive understanding of the above theorem, let 
us finally consider the tensor-network state proposed in Ref. 31 as an 
example. The state is defined by the following tensor network 31 : 

3l =|i)(o|®(i|, 



3 
~ 2 

1 

71 



= |0}{l|®<0|, 



(Z®(1| + |1)<0|®(0|), 



v/3 



(Z®(0|-|0}(1|®<1|), 
= |0}{1|, 

-|i)(o|, 



non-zero two-point correlation between a and b 



f U unitary operation 



-At-Al 



B 



a b 

no two-point correlation between a and b 

Figure 2 | Left: The AKLT state cannot be decomposed into two chains by a local unitary. Right: The tensor network of Ref. 31 . 
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B 



— |1><X| 



■|o)<o| 



(e~ 2isco cos d -A[2]+e-^ +s ^ sin^A[l]^ j ^ 
is implemented in the correlation space, where 



J cos-A[2]+e-^ +sc °) sin-A[l] 



By the assumption Eq. (3), this 



(see Fig. 2 right). A _ L 's are defined in the same way. Each horizontal 
line works as a single-qubit wire. Two nearestneighbour horizontal 
chains are decoupled by measuring sites B in the z-basis. Before 
starting the computation, the filtering operation {F, F'}, where 



F = 



F' = 



V3 



V3 



i i 



and 



2 2 

should work as a TP operation in the correlation space. Therefore, 

6 6 
yl = cos 2 -A 1 " [2] A [2]+ sin 2 -A^l^l] sinfl 

( e -^-H A t[2]A[l] + e^-^) A t[ 1 ] A [ 2 ]^ 

By the assumption Eq. (2), A 1 " [2] A [2] = £J, where £ = | |A[2] 1 1 2 . Furthermore, as we 
have shown, A f [l]A[l] = rjl (Eq. (7)). Therefore, 

y / / = e -^-5 CO ) A t[ 2 ]A[l] + e^-^)A 1 '[l]A[2], (8) 



, is applied on each site A. Let where / = \ y — £ cos 2 rj sin 2 - ] . Finally, let us consider the measurement 



us assume that the filtering is succeeded (i.e., F is realized) and a site B 

3\ 



is projected onto 



. Then, the measurement 



■{*( 



implements ZXe iZ(l) ' 2 , Xe iZ(l>/2 , or Z, respectively. (The result 

not occur.) Let us assume that the singlequdit error 
p^pEpE^ + (l-p)p, 



does 



where E is the error which exchanges 



and 



, occurs, and 



that thus affected particle is measured in M n / 2 ,^ If the measurement 
result is — ^ , we obtain the map 

p^p\l)(0\p\0)(l\ + (l-p)ZpZ, 
which is not linear CPTP, is implemented in the correlation space. 

Methods 

Proof: In order to show Theorem, let us assume that 

There is no such S. (3) 

We will see that this assumption leads to the contradiction that d < 2. First, let us 
consider the state 



, ®l/,„ 2 )| l P(L,R))f, 



(4) 



where U a ^ b = \a)(b\ + \b)(a\ +1 - \a)(a\ - \b)(b\ is the unitary error which exchanges 
\a) and | b), and I is the identity operator on a single qudit. In Eq. (4), the error U x ^ 2 is 
applied on the first physical qudit of (L,l?)) 1 . If we do the measurement Mq^ on 
the first physical qudit of Eq. (4), and if the measurement result is |2>, Eq. (4) becomes 
i d-i d-i 

£...£>|Afo,]... 

(5) 



v//»-i(|i>,A[l]|R» 



k 2 =0 k N =Q 



A{k 2 ]^^\R)\k N ,...,k 2 )®\2). 



In other words, the operation \R) - 



A[l] 



\R) is implemented in the correlation 



space. By the assumption Eq. (3), this operation should work as a TP operation in the 
correlation space. Therefore, 

At[l] A[l] _ 7 



By takings ||A[1]|| 2 , 



A t [l]A[l] = )y/. 



(6) 
(7) 



Second, let us consider the measurement Merf on the first physical qudit of 
(/® N " 1 (g)L/ 0 ^ 2 y s )|^(L,^)>f , where se{0,l, . . . ,d-l},V = J2pZ l 0 e~ iojp \p)(p\ is a 
unitary phase error, and co = 2%ld. If the measurement result is |a^), 



Mo,* on the first physical qudit of (I® N ~ l (^Uo^Uq^V'^^R)) 1 ; , where f e {0, 
1, ... y d- 1}. If the measurement result is |a^>, 

^- ! ' to cos^A[l]+e-^- 2! ' to sin^A[2]^ j Vd 



is implemented in the correlation space, where 



Vs. 



J cos-A[l] + e-^- 2ito sin^A[2] 



By the assumption Eq. (3), this should 



also work as a TP operation in the correlation space. Therefore, 

^/ = e^ +to )A t [2]A[l]+e-^ +to U t [l]A[2], 



(9) 



where S' = I S - £ sin 2 - - rj cos 2 - . From Eqs. (8) and (9), el = [ e - 2 ^-^ - 

sin 6 \ 2 2) 

gz^+to)] A t [ 2 ] a[ 1] , where e = er^y - e ^ +tco) §'. Let us assume that e~ 2 ^- s( ° > - e 2i{ * + tw) 

* 0. Then, e'l = At[2]A[l], where e* = e _^_J_ ^ +to) . If e' = 0, At[2]A[l] = 

0, which means A[l] =0 since A[2] is unitary up to a constant (assumption Eq. (2)). 
Therefore, e' ^ 0. In this case, A [1] = e"A[2] for certain e" ^ 0, since A [2] is unitary 
up to a constant. Hence e 2i ^- s(0) - e 2i{ * +tco) = 0. This means 

2(j) + (t-s)(D = r S!t n, (10) 

71 

where r s>t g {0, 1, 2, 3, . . .}. Let us take t = s = 0. Then, Eq. (10) gives (j) = r 0 ,o - (to,o e 

{0,1,2, ...}). Let us takes = 1, t = 0. Then, Eq. (10) gives = +r lQ -{r lQ e {0,1,2, 

d 2 

...}). In order to satisfy these two equations at the same time, there must exist r 0j0 and 

71 71 71 

r x o such that r 0 o - = + r\ o — ■ If r 0 0 = r x 0 , then 0 = lid which means d = °°. 

7 2 d '2 

2 

Therefore r 0 , 0 ^ r 1>0 . Then we have d = < 2, which is the contradiction. ■ 

^o,o — n,o 

One might think that if we rewrite the post-measurement state Eq. (5) as 

1 -rWMmw 



v '/n-,(|i)A[i]|«)) 1 

g.-.g(L|A[fc N ]...AN 

fc 2 =0 fc N =0 



A[l] 



r|R)|*N,...,*2>®|2) 



||A[lp)||' 

and redefine the operation implemented in the correlation space as 
A[l] 



l|A[l]|R>|| 



\R), the TP-ness is recovered in the correlation space. However, in 



this case, the non-lineally appears unless A f [l]A[l] oc J, and therefore if we require 
the linearity in the correlation space, we obtain the same contradiction. 

In short, if d > 3 not all physical errors on physical qudits appear as linear CPTP 
errors in the correlation space of pure matrix product states. 
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